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INTRODUCTION 


Several  reviewers  (Omholt,  1971,  Pemberton  and  Shepherd,  1975, 
Royrvik  and  Davis,  1977)  have  described  the  observational  character- 
istics of  quasipe  riodic  auroral  brightness  fluctuations  in  detail.  These 
authors  note  that  the  frequencies  at  which  aurorae  are  modulated  gen- 
erally fall  into  one  or  more  distinct  frequency  bands.  Fluctuations  at 
a few  Hz,  sometimes  referred  to  as  flickering,  are  often  superposed 
on  lower  frequency  pulsations  with  periods  <—  10s.  The  behavior  of 
these  two  components  is  independent,  indicating  that  the  responsible 
mechanisms  may  be  different.  The  fact  that  the  fast  pulsations  have 
periods  similar  to  the  bounce  periods  of  particles  trapped  in  the  high 
latitude  geomagnetic  field  suggests  that  the  causative  processes  in 
this  case  are  magnetospheric  in  nature.  The  present  report  is  con- 
cerned with  the  pulsations  with  periods  longer  than  the  particle  bounce 
periods. 

A model  for  the  pulsating  aurora  must  explain  all  of  the  follow- 
ing observed  features: 

1)  quasipe  riodic  brightness  fluctuations  having  pulse  shapes 
that  range  from  sinusoidal  to  half- sinusoids  which  abruptly 
"switch-on"  and  "switch-off  " 

2)  damped  pulsations;  irregular  pulsations 

3)  Pulse  amplitudes  of  up  to  100% 

4)  Pulse  periods  — 5 - 20  s (assuming  that  the  frequencies 

> 1 Hz  are  caused  by  a separate  phenomenon  as  discussed 
above) 


5)  Synchronous  pulsation  over  entire  forms  of  both  discrete 
and  diffuse  types 

6)  Streaming,  or  localized  bright  areas  that  travel  across  the 
auroral  forms. 

7)  Forms  exhibiting  an  abrupt  onset  or  cessation  of  pulsations 

8)  Preferential  occurrence  of  pulsations  at  the  e’quatorward 
edge  of  the  auroral  oval  on  closed  field  lines 

9)  Conjugate  pulsations 

1C  Preferential  occurrence  in  the  post  midnight  and  morning 
oval 

11)  A precipitating  electron  spectrum  that  varies  over  the  cycle 
of  a pulsation. 

The  theory  of  pulsating  aurorae  that  has  gained  support  in  recent 
years  is  based  on  the  hypothesis  that  pitch  angle  scattering  of  trapped 
particles  into  the  bounce  loss  cone  by  whistler  noise  is  modulated  at 
high  altitudes  by  geomagnetic  micropulsations  (Coroniti  and  Kennel, 
1970).  Recently,  the  ability  of  this  theory  to  describe  some  of  the  pro- 
perties of  the  pulsating  aurora  has  been  questioned  by  experimenters 
(Whalen  et  al,  1971.  Heacock  and  Hunsucker,  1977).  In  particular, 
Heacock  and  Hunsucker  found  that  micropulsations  appeared  to  be  the 
effect  rather  than  the  cause  of  precipitation  pulsations.  A local  or 
atmospheric  mechanism  related  to  the  ionization  caused  by  the  particles 
was  proposed  by  Maehlum  and  O'Brien  (1968)  but  they  found  that  ex- 
ceptionally intense  precipitation  was  required  to  sustain  the  self-gen- 
erated pulsations.  Observed  pulsations,  on  the  other  hand,  are  usually 
associated  with  low  intensity  aurorae  (Royrvik  and  Davis,  1977). 
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In  this  report  an  alternative  interpretation  is  presented  based  on 
observations  of  wave-like  density  fluctuations  in  the  neutral  atmosphere. 
It  is  shown  how  variations  of  the  precipitated  flux  can  be  caused  by  vari- 
ations of  the  altitude  at  which  the  particles  are  lost,  and  how  the  magni- 
tude of  the  flux  variations  is  related  to  the  form  of  the  pitch  angle  dis- 
tribution in  the  loss  cone.  The  loss  cone  variations  necessary  to  ac- 
count for  the  observed  features  of  auroral  pulsations  can  be  supplied 
by  atmospheric  acoustic-gravity  waves  similar  to  the  waves  that  dis- 
tort noctilucent  cloud  layers. 
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PRECIPITATION  MODULATIONS  FROM 
LOSS  CONE  VARIATIONS 


The  size  of  the  loss  cone  is  determined  by  the  altitude  at  which 
particles  are  lost  from  the  trapping  region  by  atmospheric  absorption 
(cf  Roederer  1970).  The  maximum  pitch  angle  in  the  loss  cone  is  de- 
fined as  the  loss  cone  angle  The  loss  cone  angle  effectively  divides 

trapped  and  precipitating  particles  in  the  magnetosphere.  From  obser- 
vations of  radiation  belt  pitch  angle  distributions  it  is  found  that  the 
pitch  angle  that  mirrors  at  ~ 100  km  altitude  typically  marks  the  bound- 
ary of  trapped  pitch  angles;  thus  reference  is  frequently  made  to  the 
"100  km  loss  cone.  " 

Consider  the  situation  illustrated  in  Fig.  1 whe.e  the  altitude  at 
which  absorption  occurs  changes  as  a result  of  some  modification  of 
the  atmosphere.  Because  the  particle  population  throughout  the  mag- 
netospheric  flux  tube  does  not  instantly  respond  to  a change  in  the  loss 
cone  at  low  altitudes,  the  distinction  between  trapped  and  precipitated 
pitch  angles  is  temporarily  upset.  In  particular,  if  the  loss  cone  in- 
creases, some  particles  that  were  previously  trapped  will  become 
members  of  the  loss  cone  population.  If  the  loss  altitude  is  steadily 
increased  over  a period  of  time,  sustained  enhanced  precipitation  will 
occur  throughout  that  period.  A similar  mechanism  is  in  operation  in 
the  vicinity  of  the  South  Atlantic  Anomaly  (Paulikas,  1975,  Torr  et  al 
1976)  where  quasitrapped  particles  drift  through  a region  of  increas- 
ing local  loss  cone  until  they  are  absorbed.  In  the  present  case,  how- 
ever, the  precipitation  occurs  by  virtue  of  a temporally  changing  local 
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Fig.  1.  Precipitation  Enhanced  by  Increasing  Loss 
Altitude  That  Produces  an  Increase  in  the 
Loss  Cone  Angle  a 


I 
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atmosphere  density  profile  rather  than  by  a spatially  changing  atmo- 
spheric density  distribution  that  is  experienced  by  azimuthally  drift- 
ing particles. 

Decreases  in  the  size  of  the  loss  cone  similarly  cause  decreases 
in  precipitation.  However,  if  the  loss  cone  increases  and  then  decreases 
sequentially,  a more  rigorous  description  of  the  magnetospheric  particle 
response  is  required  to  obtain  the  time  variations  of  the  precipitated 
flux.  Following  Torr  et.  al.  (1976),  an  equation  that  describes  the  tem- 
poral behavior  of  the  loss  cone  population  N of  a magnetospheric  flux 
tube  can  be  written 

dN  Q - L (1) 

dt 

where  Q represents  sources  and  L represents  the  loss  or  precipitation 
rate.  The  source  of  particles,  which  is  not  understood  in  detail,  pro-  « 

bably  operates  at  high  altitudes  where  the  atmosphere  has  little  in- 
fluence. If  it  is  assumed  that  all  changes  in  N are  caused  only  by  tem- 
poral changes  in  the  atmospheric  loss  cone  angle  L is  given  by 


L 


Q-(— ) 

w \ da  /«  = 


(2) 


Equation  (2)  provides  a means  of  evaluating  L,  the  rate  at  which  particles 
in  the  flux  tube  are  lost  to  the  atmosphere. 
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In  order  to  carry  out  a numerical  investigation  of  precipitation 

variations,  assumptions  regarding  the  nature  of  the  source  Q,  the 

pitch  angle  distribution  dN  and  the  rate  of  change  of  the  loss  cone 

da 

angle  da^  are  required.  According  to  Kennel  and  Petschek  (1966)  the 

nr 

introduction  of  a midplane  (ie  equatorial)  source  in  a magnetic  trapping 
geometry  wherein  pitch  angle  diffusion  is  occurring  leads  to  a mono- 
energetic  pitch  angle  distribution  at  the  loss  cone  angle  given  by 


F (*  > 


/DT)/cc 


!& T) 


where  I and  I,  are  the  modified  Bessel  functions,  S is  related  to  the 
o • 

source  strength,  T is  approximately  equal  to  a quarter  bounce  period, 
and  D is  a constant  that  characterizes  the  strength  of  diffusion.  Be- 
cause the  distribution  (3)  was  derived  under  the  assumption  of  a steady 
state,  the  quantity  S is  actually  the  rate  at  which  particles  enter  the 
loss  cone  via  the  diffusion  mechanism.  Thus,  Q in  eq.  (2)  can  be  set 

equal  to  S.  The  appropriate  value  for  /dN\  is  then  given  by  eq.  (3). 

' d a ' a = a 
c 

Note  that  S enters  eq.  (2)  as  a normalization  factor.  The  remaining 
quantity  to  be  defined  is  da 

dt~~ 
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Changes  in  the  loss  cone  size  are  caused  by  changes  in  the  al- 
titude at  which  particles  are  lost  to  the  atmosphere.  The  variations 
in  the  loss  cone  angle  related  to  a particular  perturbation  of  the  loss 
altitude  h can  be  approximated  by 


*c  (*> 


. -1 
sin 


. 2 
sin 


(4) 


where  r and  a are  values  of  the  geocentric  radius  of  the  loss  level 
o o 

and  the  loss  cone  angle  at  some  reference  point.  The  corresponding 
expression  describing  the  loss  cone  rate  of  change  da^  / dt  is 


do. 


dt 


3 

2 


tan  a 
c 


d h(t) 
dt 


l ( rQ  + h (t)) 


(5) 


In  the  present  calculations  it  is  assumed  that  rQ  corresponds  to  the 

geocentric  distance  at  100  km  altitude  and  that  aQ  has  a typical  value 

of  1.5°  which  is  the  100  km  loss  cone  at  L 9.  Three  characteristic 

loss  cone  variations  of  interest  are  : 1)  monotonically  increasing  a- 

c, 

2)  monotonically  decreasing  a and  5)  periodically  varying  a 

C t c. 

Figure  2a  shows  how  L/S  varies  with  time  for  a constantly  in- 
creasing loss  level  described  by  (6) 

h (t)  = At 

for  selected  values  of  the  parameters  D,  T and  A. 
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Fig-  2.  (a)  Variation  of  L/S  Caused  by  a Loss  Altitude  That  Constantly 

Increases  at  the  Rate  A.  The  parameter  DT  characterizes  the 
shape  of  the  pitch  angle  distribution  of  the  precipitating  parti- 
cles. (b)  Variation  of  L/S  caused  by  a loss  altitude  that  con- 
stantly decreases  at  the  rate  A. 
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It  is  seen  that  for  a given  value  of  A strong  diffusion  (large  DT)  pro- 
duces a slower  increase  in  the  loss  rate  than  weak  diffusion  (small 
DT).  This  behavior  can  be  understood  by  considering  the  shape  of  the 
loss  cone  distribution  (3).  In  strong  diffusion  the  pitch  angle  distribu- 
tion is  nearly  isotropic.  As  a result,  a specific  increase  in  the  loss 
cone  angle  will  increase  the  loss  cone  population  by  a smaller  factor 
than  the  same  increase  in  a distribution  that  rises  rapidly  outside  of 
the  loss  cone.  Pitch  angle  distributions  associated  with  weak  diffusion 
(small  DT)  possess  the  later  characteristic  (cf  Theodoridis  and  Paolini, 
1967). 

Figure  2b  shows  the  analogous  results  for  a decreasing  loss  cone 
corresponding  to  negative  A values  in  eq.  (6).  Again,  variations  are 
enhanced  when  the  value  of  DT  is  small.  In  this  case,  however,  the 
decreasing  bounce  loss  cone  can  produce  changes  in  N that  exceed  the 
rate  of  supply  by  the  source  S,  i.  e.  negative  values  of  L are  obtained 
from  eq.  (2).  Since  negative  L/S  simply  means  an  absense  of  precipi- 
tation, these  negative  values  are  set  equal  to  zero  in  Fig.  2b. 

Periodic  changes  in  the  loss  cone  size  can  be  caused  by  periodic 
changes  in  the  altitude  at  which  particles  are  lost  to  the  atmosphere. 
Under  this  condition,  variations  in  the  loss  cone  angle  can  be  obtained 
from  eq.  (4)  with 

h(t)  = A sin  nit  (7) 
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Here  A describes  the  amplitude  of  the  altitude  variation  and  cois  the 
frequency.  Equation  (5)  gives  the  corresponding  expression  for  do  / 
dt: 


do 

c 

dt 


Aid  tan  o cos  u)t/(r  + A sin  out) 
c o 


(8) 


Figure  3 shows  some  examples  of  the  periodic  variation  of  L with  time 

S 


for  selected  values  of  the  parameters  D,  T,  A and  ou.  Two  distinctive 
types  of  variations  are  apparent:  quasisinusoidal  and  pulsed.  Similar 
results  are  obtained  for  other  combinations  of  D,T,A  and  ou.  These 
two  forms  can  be  qualitatively  explained  by  considering  the  dynamic 
competition  between  high  altitude  diffusion  into  the  loss  cone  and  the 
"windshield  wiper"  action  of  the  atmosphere  at  the  ends  of  the  flux  tube. 
If  practically  sinusoidal  modulation  occurs,  the  diffusion  must  be  strong 
enough  to  produce  a nearly  isotropic  pitch  angle  distribution.  The  vari- 
ation of  the  loss  cone  population  then  virtually  tracks  the  loss  cone  os- 
cillations at  low  altitude.  If  pulsed  modulations  occur,  diffusion  is 
sufficiently  weak  that  a significant  density  gradient  exists  at  the  loss 
cone  angle.  Precipitation  is  enhanced  during  the  increasing  phase  of 
the  loss  cone  cycle,  and  decreases  as  the  loss  cone  shrinks,  but  in  the 
latter  phase  the  diffusion  rate  cannot  replenish  the  newly  trapped  pitch 
angles  fast  enough  to  maintain  the  equilibrium  distribution  (3).  As  be- 
fore, when  the  diffusion  rate  into  the  loss  cone  cannot  keep  up  with  the 
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Fig.  3.  Variation  in  L/S  Produced  by  Sinusoidal  Variations  in  the  Loss 
Altitude.  These  pulsations  occur  for  a variety  of  wave  ampli- 
tudes A,  frequencies  u>,  and  diffusion  parameters  DT.  Nega- 
tive values  of  L/S  have  been  set  equal  to  zero. 


rate  at  which  loss  cone  particles  can  be  lost,  equation  (2)  gives  negative 
values  of  L which  have  been  set  equal  to  zero  in  Fig  3 to  indicate  an  ab- 
sense  of  precipitation.  The  behavior  in  this  case  is  similar  to  the  response 
created  by  the  rapid  withdrawal  of  a piston  in  a container  of  gas;  the  strong 
diffusion  case  would  correspond  to  a slowly  withdrawn  piston. 

The  energy  spectrum  of  precipitating  particles  will  vary  during  ' 

the  cycle  of  a pulsation  if  the  magnitude  or  form  of  the  pulsations  depends 

on  the  energy.  Figure  4 illustrates  how  either  spectral  hardening  or  soften-  ' 

t 

ing  can  occur  at  the  peak  of  a pulsation.  Because  both  the  diffusion  coef- 

I 

ficient  D and  the  time  constant  T depend  on  tne  particle  energy,  differences  1 

i 

in  monoenergetic  pulse  shapes  and  hence  spectral  variations  are  an  ex-  1 

pected  consequence  of  the  present  model.  J 

In  reality,  the  source  S and  pitch  angle  distribution  dN  may  not  \ 

da 

maintain  their  equilibrium  values  (cf  eq.  (3))  in  all  of  the  cases  under  1 

consideration.  However,  an  exact  treatment  would  require  the  non- 

« 

trivial  solution  of  the  time  dependent  diffusion  equation  with  a time 
dependent  boundary  condition  at  a = a^.  The  qualitative  behavior  shown 
in  Figs  2-3  is  probably  a good  approximation  provided  that  no  instabilities 
associated  with  the  loss  cone  variations  modify  the  diffusion  coefficient 
or  the  source. 

ATMOSPHERIC  WAVES  AS 
SOURCE  OF  LOSS  CONE  VARIATIONS 

| 

I 

Atmospheric  pressure  waves  have  been  observed  since  the  early 
sixties  (cf  Chrzanowski  et  al  1961,  Campbell  and  Young,  1963).  Wilson 
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E1  h E3 


Fig.  4.  (a)  Three  Different  Precipitation  Pulse  Shapes 

at  Three  Energies,  (b)  Particle  spectrum  be- 
fore pulse,  at  pulse  peak,  and  after  pulse  when 
the  energies  Ej»  E^,  and  E3  are  ordered  as 
shown.  The  spectrum  softens  at  the  pulse  peak, 
(c)  Particle  spectrum  before  pulse,  at  pulse 
peak,  and  after  pulse  when  energies  Ej,  E^> 
and  E3  are  ordered  as  shown.  The  spectrum 
hardens  at  the  pulse  peak. 
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(1967)  proposed  a shock  wave  model  for  the  generation  of  infrasonic 
waves  by  supersonically  moving  auroral  forms.  At  about  the  same 
time,  the  existence  of  longer  period  waves  in  the  neutral  atmosphere 
was  established  from  observations  of  traveling  ionospheric  disturbances 
(TIDs)  (eg.  King,  1967,  Bowman,  1968).  Within  the  last  decade , theo- 
retical models  have  linked  the  source  of  these  neutral  waves  with  the 
aurora  (Thome,  1968,  Davis  and  da  Rosa,  1969,  Chimonas  and  Hines, 

1970).  General  theories  of  waves  in  the  neutral  atmosphere,  the  acoustic- 
gravity  waves,  have  since  received  considerable  attention  (Hines,  1974, 
Francis,  1975,  Smith,  1977). 

Atmospheric  waves  at  frequencies  above  and  below  ~ 10  ^s  * be- 
long to  the  acoustic  and  gravity  domains,  respectively.  The  higher  fre- 
quency acoustic  waves  are  generally  attenuated  within  a few  100  km  of 
the  source,  while  gravity  waves  can  propagate  distances  the  order  of  an 
earth  radius  (Georges,  1968).  However,  acoustic  waves  also  appear  to 
propagate  in  atmospheric  "sound  channels"  to  mid-latitudes  from  the 
auroral  zone  (Chrzanowski  et  al.  , 1961).  The  picture  of  neutral  waves 
that  emerges  from  integrating  theory  and  observation  (Newton  et  al.  , 

1969,  Reber  et  al.  , 1975)  is  one  in  which  wave-like  corrugations  of  the 
isodensity  surfaces  of  the  atmosphere  above  ~ 80  km  are  generated  by 
active  auroral  arcs  and  launched  approximately  along  magnetic  meridians 
with  horizontal  phase  velocities  equal  to  the  local  sound  speed  of  ~ 200  m 
s'1.  The  wave  amplitude  expressed  in  terms  of  departure  from  the 
quiescent  atmospheric  density  as  given  by  standard  model  atmospheres 
is  typically  — 15%  - 50%.  Noctilucent  cloud  observations  (Fogle  and  Haurwitz, 
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1966)  indicate  the  presence  of  horizontal  wavelengths  at  ~ 85  km  al- 
titude of  ~ 1 - 10  km,  while  horizontal  wavelengths  of  ~ 100  - 1000  km 
are  inferred  from  TID's  (Francis,  1975).  The  vertical  extent  or  am- 
plitude of  these  waves  covers  a wide  range,  but  is  often  comparable  to 
the  horizontal  wavelength. 

Although  the  periods  of  gravity  waves  are  much  too  long  to  ex- 
plain the  frequencies  of  typical  auroral  pulsations,  acoustic  waves  with 
periods  ~ 10s  are  also  present  in  the  auroral  zone  atmosphere.  The 
~ 1 - 10  km  wavelengths  that  are  seen  traveling  horizontally  at  the  sound 
speed  in  noctilucent  cloud  layers  can  provide  the  observed  auroral  pul- 
sation frequencies.  Moreover,  as  illustrated  by  Fig.  5 . thermospheric 
winds  which  blow  across  the  polar  cap  at  velocities  in  excess  of  ~ 100  m 
s ^ (Meriwether  et  al.  , 1973,  Brekke  et  al.  , 1974)  can  Doppler  shift 
neutral  waves  to  frequencies  between  ~ 0 and  twice  their  natural  fre- 
quency. 

DISCUSSION 

In  view  of  the  thesis  of  the  paper,  it  is  interesting  to  note  that 
Porter  et  al.  (1974)  pointed  out  the  possible  connection  between  neu- 
tral waves  and  airglow  intensity  fluctuations.  Auroral  pulsations,  on 
the  other  hand,  have  previously  been  considered  as  a propective  source 
of  neutral  waves  (Maeda  and  Watanabe,  1964). 

In  the  introduction  a list  of  the  observational  properties  of  the 
pulsating  aurora  was  given.  This  set  of  characteristics  can  now  be 
reconsidered  from  the  viewpoint  of  the  present  model.  In  fact,  all  of 
the  listed  observational  features  can  be  explained: 
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Fig.  5.  Doppler  Shift  of  Neutral  Waves  by 
Polar  Cap  Neutral  Winds.  The 
phase  and  wind  velocities  cancel  on 
the  dayside,  while  the  velocities 
reinforce  on  the  nightside  where  pul- 
sations are  observed. 
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(1)  Figure  3 illustrates  that  both  sinusoidal  and  pulsed  pre- 

cipitation can  be  driven  by  the  passage  of  a sinusoidal 
wavetrain. 

(2)  Damped  pulsations  correspond  to  damped  waves;  irregular 
pulsations  can  result  from  aperiodic  local  density  fluctuations. 

(3)  Pulse  amplitudes  within  the  framework  of  this  model  are 
virtually  unlimited. 

(4)  The  periods  of  the  precipitation  variations  are  limited  only 
by  the  frequency  spectrum  of  the  atmospheric  density  fluctu- 
ations and  by  the  magnitude  and  direction  of  the  local  neutral 
wind  velocity. 

(5)  An  entire  form  can  pulsate  synchronously  if  the  wavefront 

lies  parallel  to  the  long  axis  of  the  form  as  illustrated  in  Fig.  6a. 

(6)  Figure  6b  shows  how  streaming  can  occur  if  the  wave  pattern 
propagates  at  a large  angle  to  the  long  axis  of  the  form. 

(7)  Pulsations  can  turn  on  and  off  without  disturbing  the  mor- 
phology of  the  form  as  the  finite  length  wavetrain  passes. 

(8)  Waves  generated  at  high  latitude  where  active  arcs  are 
found  will  cause  pulsations  of  lower  latitude  arcs  and  diffuse 
forms. 

(9)  Pulsations  will  occur  conjugately  (within  a bounce  period) 
if  backscatter  is  sufficiently  intense  and  if  the  conjugate 
point  atmosphere  is  not  undergoing  comparable  oscillations. 

(10)  The  polar  cap  neutral  wind  pattern  usually  shows  a strong 

afternoon  quadrant  to  predawn  quadrant  alignment  (cf  MeriT 
wether  et  al,  1973);  only  waves  traveling  equatorward  in 
the  vicinity  of  the  predawn  quadrant  will  be  Doppler  shifted 
to  the  frequencies  that  characterize  auroral  pulsations. 

(11)  Precipitation  caused  by  the  oscillating  atmosphere  depends 
on  the  particle  diffusion  coefficient  D and  quarter  bounce 
time  T (ie.  on  the  pitch  angle  distribution  determined  at 

high  altitude)  which  are  energy  dependent;  thus,  the  pulsations 
will  be  different  at  different  energies.  The  largest  fractional 
fluctuations  are  associated  with  weak  diffusion  or  mirroring 
pitch  angle  distributions.  Distributions  peaked  at  90  have 
been  observed  by  Whalen  et  al  (1975)  near  a pulsating  aurora. 
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SYNCHRONOUS  STREAMING 

PULSATION 


Fig.  6.  (a)  Synchronous  Pulsations  That  Can  Occur  If  the  Wavefront 

Lies  Parallel  to  the  Auroral  Arc.  (b)  Streaming  that  can  be 
caused  by  a wavefront  traveling  perpendicular  to  the  long 
axis  of  an  auroral  form. 
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The  model  for  the  pulsating  aurora  proposed  here  is  by  no  means  com- 
plete. In  particular,  the  morphology  of  neutral  waves  at  ~ 85  - 100  km 
altitude  where  the  auroral  emissions  are  generated  can  only  be  inferred 
from  observations  of  noctilucent  clouds.  However,  as  indicated  by  the 
foregoing  list,  the  model  can  at  least  conceptually  account  for  a large 
number  of  the  experimentally  determined  characteristics  of  auroral  pul- 
sations. 
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THE  IVAN  A.  GETTING  LABORATORIES 

The  Laboratory  Operation*  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  cnemi cal  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro- optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix_con5osTtes”an3^iewToritis~of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magneto  spheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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